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Aims. By combining the VIMOS VLT Deep Survey (VVDS) with the Spitzer Wide-area InfraRed Extragalactic survey (SWIRE) data, we have 
built the currently largest spectroscopic sample of high redshift galaxies selected in the rest-frame near-infrared. In particular, we have obtained 
I . 2040 spectroscopic redshifts for a sample of galaxies with a magnitude measured at 3.6/um (m 3 6 ) AB < 21.5, and 1255 spectroscopic redshifts for 
a sample of galaxies with (ni4 5 ) AB < 21. These allow us to investigate, for the first time using spectroscopic redshifts, the clustering evolution of 
£ — ' galaxies selected from their rest-frame near-infrared luminosity, in the redshift range 0.2 < z < 2.1. 

f"^ , Methods. We use the projected two-point correlation function w p (r p ) to study the three dimensional clustering properties of galaxies detected 
— ■ at 3.6/um and 4.5pm with the InfraRed Array Camera (IRAC) in the SWIRE survey with measured spectroscopic redshifts from the first epoch 
• VVDS. In addition, we measure the clustering properties of a larger sample of 16672 SWIRE galaxies for which we have accurate photometric 
Or redshifts on the same field by computing the angular correlation function. We compare these measurements. 

Results. We find that in the flux limited samples at 3.6pm and 4.5yum, the apparent correlation length does not change from redshift ~ 2 to the 
present. The measured correlation lengths have a mean value of r = 3.9 ± 0.5 h~ [ Mpc for the galaxies selected at 3.6/jm and a mean value of 
ro =! 4.4 ± 0.5 h~ l Mpc for the galaxies selected at 4.5/im, all across the redshift range explored. These values are larger than those typicaly found 
^ ' for I-band selected galaxies at I AB < 24, for which r varies from 2.69 \r x Mpc to 3.63 /T 1 Mpc between z = 0.5 to z = 2.1. We find that the 
. . , difference in correlation length between I-band and 3.6 - 4.5/um selected samples is decreasing with increasing redshift to become comparable 
^ . at z — 1.5. We interpret this as evidence that galaxies with older stellar populations and galaxies actively forming stars reside in comparably 
over-dense environments at epochs earlier than z — 1.5, supporting the recently reported flattening of the color-density relation at high redshift. 
The increasing difference with cosmic time in correlation length observed between rest-frame UV-optical and near-infrared selected samples could 
JH , then be an indication that star formation is gradually shifting to lower density regions as cosmic time increases, while the older passively evolving 
. . i galaxies remain to trace the location of the highest primordial peaks. 

Key words. Cosmology: observations - Cosmology: large scale structure of Universe - Galaxies: evolution - Galaxies: high-redshift - Galaxies: 
statistics - Infrared: galaxies 
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1. Introduction 

According to the current cosmological paradigm, the formation 
of large scale structures can be described as the evolution of pri- 
mordial dark matter mass density perturbations under the influ- 
ence of gravity. The galaxies that we observe are the result of 
the cooling and fragmentation of gas within the potential wells 
provided by the dar k matter halos, that hierarchically build-up 
(IWhite & ReeslH978l) . Even if the distribution of galaxies is in 
some way biased w ith respect to the distribution of the dark 
matter (Kaiser 1984), the spatial disnibution of galaxies should 
therefore trace the dark matter density field. It is to be expected 
that the physical processes that build galaxies of different types 
and luminosities are sensitive to the mass of halos and their dif- 
ferent environments. Therefore the evolution of the clustering 
of galaxies may depend on galaxy type, luminosity and envi- 
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ronment. The comparison of the clustering of different galaxy 
populations along cosmic time could therefore provide useful 
constraints on the evolution and formation scenario of galaxies. 
When galaxies are observed from redder wavelengths and up to 
a few microns rest-frame, this preferentially selects older stellar 
populations formed earlier in the life of the Universe. As the old- 
est stars are likely to have formed in the highest density peaks in 
the Universe, it is therefore expected that the clustering of their 
host galaxies would be significantly higher than that of galaxies 
hosting more recently formed stars. In the local Universe, se- 
lecting galaxies at 2. 15 micr ons (K s band) fr om the Two Micron 
All Sky Survey (2MASS, iMaller et al.ll2005h . it is found that 
these galaxies are more clustered than those selected at opti- 
cal wavelengths, with an angular correlation amplitude several 
times larger than found in t he Sloan Digital Sky Survey Early 
Data Release (SDSS EDR, IConnollv et al J 120021) , even if part 
of this difference may be due t o the different lumi nosities of 
the samples. At higher redshifts, Oli ver et a fl (12004 . based on 
Spitzer data, found that at Zmedkm — 0-75 the correlation length of 
(m3.6)AB < 20.1 (or equivalently Kv ega < 18.7) selected galaxies 
is ro = 4.4 + 0.1 h~ l Mpc with an angular correlation amplitude 
~ 60 times lower than for 2MASS local galaxies, and the more 
lumin ous galaxies are even more strongly clustered (Farra h~et alj 
120061) . 

The IRAC 3.6pm and 4.5pm bands are only marginally affected 
by the first polycyclic aromatic hydrocarbon (PAH) spectral fea- 
ture at 3.3pm, which fall in these bands in the redshift ranges 
[0,0.18] and [0.21,0.51] respectively. However, galaxy popula- 
tions observed in these bands contain a mix of early-type, dust- 
redde ned, and star-forming systems dRowan-Robinson et alj 
120051) which make clustering expectations less straightforward. 
This paper aims to characterize the evolution of the real-space 
correlation length ro of near-infrared selected galaxies from 
the Spitzer Wide-are a InfraRed Extragalactic survey (SWIRE, 
lLonsdale et al.120031) . which have secure spectro scopic redshifts 
from the VIMOS VLT Deep Survey (VVDS, ILe Fevre et all 
2005b]). We compare these mea surements to those in shallower 
samples, e.g. lOliver et alj d2004l) . as well as to th e observed evo- 
lution o f optically-selected galaxies as shown in lLe Fevre et alj 
d2005al) . up to z - 2. 

In Section 2 we describe the galaxy sample that we use for this 
analysis. We present in Section 3 the formalism of the two-point 
correlation function and the results. Finally in Section 4 we dis- 
cuss the observed clustering evolution and we conclude. 
Throughout this analysis we assume a flat ACDM cosmology 
with Q M = 0.3, Q A = 0.7, cr 8 = 0.9 and H = 100 h km ■ s~ x ■ 
Mpc^ 1 , while a value of Hq = 70 km ■ s • Mpc~ l is used when 
computing absolute magnitudes. 

2. The VVDS-SWIRE sample 

The VVDS-SWIRE field is th e intersection of the VIRMOS 
Deep Imaging Survey (VDIS, iLe Fevre et al.1 120041) and the 
SWIRE-XMM/LSS fields, the total area covered on the sky is 
~ 0.82 deg 2 . From this field we extract two galaxy samples. 
The largest one consists of those SWIRE galaxies present in 
the field for which we measured photometric redshifts, called 
in the following the photometric redshift sample. The technique 
that we use to measure photometric redshifts is presented in 
the next section. The other sample, hereafter the spectroscopic 
redshift sample, contains galaxies included in a sub-area of 
~ 0.42 deg 2 which correspond s to a fraction of the VVDS-Deep 
survey dLe Fevre et al. 2005b) area and for which we have spec- 
troscopic redshifts measured during the VVDS first epoch obser- 



vations. We restrict this sample to galaxies with secure redshifts, 
i.e . with a redshift confid ence level greater than 80% (flag 2 to 
9. ILe Fevre et alj|200 5b). About 25% of the galaxy population 
is randomly sampled and we are correcting the correlation mea- 
surements for geometrica l effects, e.g. on s mall scales, as de- 
scribed in Section 3 and in lPollo et alJd2005l) . In these two sam- 
ples the star-galaxy separation has been performed first by using 
optical and mid-infrared colors to define the locus of galactic 
stars. In addition we used the spectral energy distribution (SED) 
template fitting procedure to reject any remaining stars and QSO. 
Objects which are best fitted with a star or a QSO SED tem plate 
and which have a high SExtractor (Ber tin & Arnouts| [l996) stel- 
larity index in the i' band (CLASS_STAR > 0.97) have been re- 
moved from the samples. The full description of the data-set and 
its properties is given in Le Fevre et al. (2007, in preparation). 

2. 1 . Photometric redshifts 

The photometric redshifts have been estimated using the x 2 fit- 
ting algorithm Le Phar$\ and calibrated with the VVDS first 
ep och spectros c opic r edshifts on the same data-set as described 
in Ilber t et all (|2006), 

including in addit ion the 3.6pm an d 
4.5pm infrared photometry from SWIRE (Lons dale et al.l l2003). 
The photometric data cover the wavelen gth domain 0.3pm < 
A < 5.0pm including B,V,R,I from VDIS £e Fevre et al.ll2004 . 
u*,g',r',i',z' from the Canada France Hawaii Telescope Legacy 
Survey (CFH TLS-D1, McCrack en et al., in preparation), K-band 
from VDIS (llovino et alj I2005T) an d the 3.6pm and 4.5p m in- 
frared photometry from S WIRE dLonsdale et alj [2003). The 
empirical templates used by lllbert et alj (l2006l) have been ex- 
tended to the mid-infrare d domain using the GISSEL library 
(Bruz ual & Charlotll2003h and we use the spectroscopic sample 
to simult aneously adapt th e templates and the zero-point cali- 
brations (lllbert et al.ll2.006l) . The comparison with 1500 spectro- 
scopic redshifts from the VVDS first epoch data shows an ac- 
curacy of cr Az / (1+Z ) 0.031 with a small fraction of catastrophic 
errors: 1.7% objects with |Az| > 0.15(l+z), up to a redshift z ~ 2. 
The full description of the method and t he data used to com pute 
these photometric redshifts is given in lllbert et all (120061) and 
Arnouts et al. (2007, in preparation). 

2.2. Selection functions 

We split the spectroscopic redshift and photometric redshift sam- 
ples into two sub-samples, applying a simple magnitude selec- 
tion from observed fluxes on the two first infrared photomet- 
ric bands of the SWIRE survey, the Spitzer-IRAC 3.6pm and 
4.5pm bands. These sub-samples satisfy respectively to a magni- 
tude in the AB system measured in the 3.6pm band of (m^ ,(,)ab < 
21.5 and in the 4.5pm band of (m^AB < 21 or, equivalently in 
flux, S3.6 > 9.2 pjy and S 4.5 > 14.5 pJy. These criteria corre- 
spond to the limits of completeness for the SWIRE bands (Le 
Fevre et al., 2007, in preparation). In addition, the spectroscopic 
redshift sample has a Iab cut due to the fact that the VVDS-Deep 
spectroscopic sample is a magnitude-selected sample according 
to the criterion 17.5 < Iab < 24. The Iab cut introduces a color 
dependent incompleteness in terms of SWIRE magnitudes in the 
spectroscopic redshift sample. Indeed we loose a small fraction 
of objects which correspond to faint and red objects as shown 
in FigureQ]and Figure[2]for the 3.6pm and 4.5pm galaxies. We 
estimate the fraction of lost galaxies to be respectively 11.6% 
and 10.3%. To understand the population missed in the spec- 

1 http://www.oamp.fr/arnouts/LE_PHARE.html 
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Fig. 1. Color-magnitude diagram for the 3.6fim spectroscopic 
and photometric sub-samples. All the galaxies detected at 
3.6fim in the VVDS-SWIRE field are represented in grey scale. 
The photometric redshift sub-sample consists of all the objects 
with (;«3.6)ab < 21.5 and galaxies from the spectroscopic sub- 
sample are shown as crosses. The vertical line indicates the limit 
of completeness and the two diagonal lines show the limits of 
the spectroscopic selection 17.5 < Iab < 24. We estimate that 
the fraction of galaxies we lose in the spectroscopic redshift sub- 
sample due to the VVDS spectroscopic selection and assuming 
the completeness limit {m^^AB < 21.5, is 11.6%. 



Fig. 2. Color-magnitude diagram for the 4.5/im spectroscopic 
and photometric sub-samples. All the galaxies detected at 
4.5^/w in the VVDS-SWIRE field are represented in grey scale. 
The photometric redshift sub-sample consists of all the objects 
with {m^AB < 21 and galaxies from the spectroscopic sub- 
sample are shown as crosses. The vertical line indicates the limit 
of completeness and the two diagonal lines show the limits of 
the spectroscopic selection 17.5 < Iab < 24. We estimate that 
the fraction of galaxies we lose in the spectroscopic redshift sub- 
sample due to the VVDS spectroscopic selection and assuming 
the completeness limit (m^AB < 21, is 10.3%. 



troscopic sample, we have used the photometric redshifts de- 
scribed in the previous section. We find that the galaxies that we 
are missing in the 4.5/j.m spectroscopic redshift sub-sample are 
redder than those in the 3.6/im spectroscopic sub-sample. The 
median (J -m-i ^AB color of these sub-samples are 3.72 and 3.67 
respectively, as shown in Figure[3] We note that photometric red- 
shifts may have difficulties in identifying dust obscured AGN 
because we are lacking the redder bands like the Spitzer-IRAC 
24/jm. Ho wever, the fraction of dust reddened AGN is expected 
to be low (iTajer et a l.ll2007i), and its impact on clustering prop- 
erties must be insignificant. Therefore, the photometric redshift 
sample, which is assumed to be complete and free of significant 
contaminants, allows us to evaluate the effect of the color depen- 
dent magnitude incompleteness on the clustering measurements 
performed on the spectroscopic redshift sample. 

3. Clustering measurements 

3.1. Methods 

3.1.1. Computing the projected correlation function 

The method applied to compute the projected correlation func- 
tion and to derive the real-space correlation length ro and the cor- 
relatio n function slope y, is described extensively in Poli o et all 
(2005). We present in this section the main points of the method. 
First we evaluate the bi-dimensional two-point correlation func- 
tion £(r p ,7t) using the standard Landy and Szalay estimator 



(lLandv & Szalaviri993l) . This estimator of the two-point corre- 
lation function is defined by: 

_ GG(r)-2GR(r) + RR(r) 

t {r) = Mr) (1) 

where GG(r), RR(r) and GR{r) are the normalized numbers of 
independent galaxy-galaxy, galaxy -random and random-random 
pairs with comoving separation between r and r + dr. The three- 
dimensional galaxy space distribution recovered from a redshift 
survey and its two-point correlation function are distorted 
due to the effect of galaxy peculiar velocities. Therefore the 
redshift-space separation s differs from the true physical comov- 
ing separation r. These distortions occur only radially since pe- 
culiar velocities affect only redshift. Thus we compute the bi- 
dimensional two-point correlation £(r p , n) by splitting the sepa- 
ration vector s into two components: r p a nd n, respectively p er- 
pendicular and parallel to the line of sight (Fisher et al. 1995- In 
this way, we separate the redshift distortions from the true spa- 
tial correlations. 

Then, we project £;(r p ,n) along the line of sight, onto the r p axis. 
This allows to integrate out the dilution produced by the redshift- 
space distortion field and we obtain the projected two-point cor- 
relation function, which is defined by, 

= 2 J o f(^rl+f}dy (2) 

In practice, the upper integration limit has to be finite to avoid 
noise, therefore we choose its optimal value as 20 h Mpc 
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Fig. 3. The right panels present the color-magnitude diagrams 
(I - m-i ^AB versus (m^ ^AB for the 3.6/um (top) and 4.5yum (bot- 
tom) photometric sub-samples considering the objects that fall 
above the spectroscopic selection upper limit Iab = 24. The 
galaxy number density is shown in grey scale and the dashed line 
indicates the median of the {I - m^AB distributions. In the left 
panels are shown the angular correlation function for these two 
selections in the redshift range z P hot = [1, 1.7]. In these panels 
the solid lines show the best-fit model when letting the two pa- 
rameters free, while the dashed lines indicate the best-fit model 
with fixed y = 1.8. 



(iPollo etalJ 12005). Assuming a power-law form for ^(r), i.e. 
£(r) = (r ( )/r) r , one can write w p (r p ) as, 



w p (r p ) = r p — 



>•(> 



(3) 



where F is the Euler Gamma function. 

We fit the w p (r p ) measurements to a power-law model (equa- 
tion 3) by minimizing the generalized^ 2 . The advantage of this 
method is that it takes into account the fact that the different 
points of the correlation func tion are correlate d. We define the 
generalized^ in our case as (Polio et al. 200J), 
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N„ N„ 



wTV,) 
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where w p (r ; ) (w° fa)) is the value of W p 



(w° p hs ) computed 

at r p = r-i with 1 < i < N p (N p is the number of r p points), 
and C is the covariance matrix of the data estimated using boot- 
strap re-sampling. 100 bootstrap samples have been used in this 
analysis. The error bars estimated from the covariance matrix 
associated to the bootstrap samples include the statistical er- 
ror only and do not take into account the error due to cosmic 
variance. Therefore, in order to have a realistic estimation of 
the errors bars on our measurements, we compute w p (r p ) on 50 
mock VVDS surveys constructed using the GallCS simulation 



dHatton et al.l l2003). This simulation of hierarchical galaxy for- 
mation is based on a hybrid N-body/semi-analytic model and is 
particularly well adapted to simulate the propertie s of the high- 
redsh i ft galaxies that w e observe in the VVDS dBlaizot et alj 
120051: IPollo et alj|2005l) . Moreover, the selection function and 
the observational biases present i n the VVDS dat aset have been 
included in these mock samples (Polio et al. 2005). Thus we fit 
the w p (rp) measurements from the simulations using the same 
method as described previously and we estimate the confidence 
limits on the ro and y measurements lo oking at the Ay 2 distribu- 
tion from the best fit (IPollo et al.ll2005l) . The error bars obtained 
using this method are larger but more realistic than those from 
the bootstrap method. They include not only the statistical error 
but also the error which comes from the cosmic variance, which 
dominates here. 



3.1.2. Computing the angular correlation function 

In order to compute the angular two-point correlation function 
w(0) (ACF), we use the Landy & Szalay estimator presented in 
equation [1] but in this case, simply considering angular pairs. 
Assuming that w{6) is well described by a power-law model as 
w(6) = A w 9 l ~ y ', we fit the measured ACF to this model using 
the generalized^ 2 method and deduce the amplitude A „. and the 
slope 1 —y. Because of the finite size of the survey, we introduce 
the integral constrain t (IC) correction in our fitting procedure. It 
is defined as follows (Roch e etaljf 19931) . 



IC 



with, 



Wy) 



e l - y dciidn 2 



(6) 



where Q. is the area of the observed field. We compute B(y) by 
numerically integrating this expression over the entire data area 
excluding regions where there are photometric defects. At the 
end our measurements are fitted with the expression: 

HW(0) = A w l ~ y - IC = A H ,(6 l - 



-y _ 



B{y)) 



(7) 



In order to derive the real-space correlation length ro from the 
angular correlation amplitude A w , we use the Limber deprojec - 
tion technique as described in Maglioc chetti & Maddoxl (fl999). 
This method, w hich is based on the Limber relativistic equation 
(Limber 1953), permits to recover the value of ro given A w and 
the redshift distribution of the galaxies. In general, this technique 
is not very stable and the shape of the redshift distribution that 
one considers is critical, but in this study, we measure ACF in 
relatively narrow redshift bins so the shape we assume for the 
inversion is not so critical. We use a smoothed redshift distribu- 
tion by applying a moving window with a width of Az = 0.2 
to the initial photometric redshift distribution (Le Fevre et al., 
2007, in preparation). The width of the smoothing window has 
been chosen to reflect the typical uncertainty that we have on 
the photometric redshift values. The ACF error bars have been 
estimated using the bootstrap re-sampling method. 

3.2. Clustering results 

3.2.1 . Clustering of the spectroscopic redshift sample 

We compute w p (r p ) in increasing spectroscopic redshift slices 
from z — 0.2 to z — 2.1. Therefore we divide the 3.6fim sub- 
sample in six redshift slices: [0.2,0.5], [0.5,0.7], [0.7,0.9], 
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Table 1. VVDS-SWIRE spectroscopic redshift sample: sub-sample properties and associated measurements of the slope, the corre- 
lation length and the linear bias. 

Spectroscopic redshift sample 



Sub-sample 


Redshift 
interval 


Mean 
redshift 


Number 
of galaxies 


(Mb) Vega 


r 


r„ (hr l Mpc) 


r (h 1 Mpc) 
(fixed y = 1.8) 




3.6pm 


[0.2, 0.5] 


0.36 


336 


-19.53 


1 92+0-27 


3 " 8+ Q7 


3 9 +u 6 

4 4+8:? 

3 7+0 + 
43+8:t 


93 +ulu 




[0.5,0.7] 
[0.7, 0.9] 
[0.9, 1.1] 


0.60 
0.81 
0.99 


445 
512 
385 


-20.29 
-20.75 
-21.17 


1 77+0.15 
j g7+8:H 

2 44+8:2! 
7 i Q +8:i9 

^ -0.38 

1 84 +0 18 
1 84 +014 

i O7+0.26 

t.»z_ 022 


4.2t| 

3.8!; 
42+0.4 


115^ 
" -Q.1Q 

1 i c+0.09 
• -0.09 

1 38 


4.5pm 


[1.1,1.3] 
[1.3,2.1] 
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[0.7, 1.1] 
[1.1,2.1] 


1.19 
1.57 
0.48 
0.90 
1.42 


195 
93 
485 
537 
182 


-21.66 
-22.13 
-20.22 
-21.15 
-21.97 


a* 

t q+0.5 
J - y -0.5 

4 +03 

t-O_ 04 
o 7+0.6 
3 - '-0.7 


t 7+8:1 
76^ 

'•"-1.0 
t q+0.5 

4 +0 3 

a 7+0.6 
J -'-0.7 


1.62^ :4 | 
7 cc+8:^i 

-096 

i or" 08 

-0.08 

1 26+ 005 
1,z - u -o.io 

1 48 +024 
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Fig. 4. ILeft panels: projected correlation functions for the 3.6pm spectroscopic redshift sub-sample in different redshift slices with 
their best-fit power law. Right panels: the 68.3%, 90% and 95.4% confidence contours for the clustering parameters ro and y. 



[0.9, 1.1], [1.1, 1.3], [1.3,2.1]. The 4.5//m sub-sample has fewer 
objects than the 3.6pm one because of the slightly shallower 
depth. In order to maximize the number of objects per redshift 
slice, we consider in this case three slices only: [0.2,0.7], 
[0.7, 1.1] and [1.1,2.1]. 

In the fitting procedure, we use all the w p (r p ) points for 
0.1 < r p < 12 h~ l Mpc. The projected correlation functions 
computed in each redshift slice are plotted in Figure [4] with the 
contours of confidence for the best-fitted parameters. Table Q] 
summarizes the values of the correlation length and the corre- 
lation function slope that we derive for the two sub-samples of 
the spectroscopic redshift sample. 



The correlation length measured in the 3.6pm and 4.5pm sub- 
samples is roughly constant over the redshift range z=[0.2,2.1]. 
It varies slightly between 3.8 and 4.2 lr l Mpc, except in the 
high-redshift interval [1.3,2.1] of the 3. 6pm sub-sample where a 
large correlation length of 5.8 h~ l Mpc is found. In this redshift 
range we are probing the intrinsically brightest galaxies, with 
a median absolute magnitude of (MB)vega = -22.13 (Table 
[TJ. Compared to the ch aracteristic absolu te magnitude M* at 
this redshift reported by lllbert et al] d2005l) . it is found that for 
these objects, AM* = (M B ) - M* =s -0.27 and equivalently 
(Lb)/L* b ^ 1.3. Thus, these luminous galaxies are expected 
to be more strongly clustered as the halo clustering strength 
increases more rapidly for luminosities greater than a few times 
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L* dZebavietd . 2005; Norberg et alj|2002l) . Indeed iPollo et al J 
(2006l): ICoil et al.l d2006) have shown that the galaxy correlation 
length increases more steeply for magnitude brighter than 
M* at z — 1 with ro values about 5 h Mpc. Our high ro 
measurement in this high redshift bin is thus fully consistent 
with these findings. Moreover, in this redshift interval the 
1 .(spin rest-frame emissi on feature from the photosp heric emis- 
sion of evolved stars dSimpson & Eisenhardtlll999l) . is present 
in the 3.6pm band, favoring the detecti on of more massiv e and 
significantly more clustered systems (Far rah et al.l 12006). We 
discuss in details the evolution of the correlation length for the 
two sub-samples in Section 4. 



3.2.2. Clustering of the photometric redshift sample 

We measure the ACF and deduce the clustering parameters ro 
and y on the photometric redshift sample in order to test the mea- 
surements from the spectroscopic redshift sample and to evaluate 
the effect of incompleteness. 

With the same methodology we compute the ACF in increas- 
ing photometric redshift slices, selecting the slice boundaries in 
order to take the photometric redshift errors into account. We 
consider five slices with increasing size: [0.2,0.5], [0.5,0.7], 
[0.7, 1.0], [1.0, 1.3] and [1.3, 1.7]. No redshift slices have been 
considered beyond z p hot = 1 .7 because the lack of accuracy of 
the photometric redshifts in this case seriously corrupts the mea- 
surements. 

We fit the measured ACF to a power-law model over 0.0017 < 
6> < 0.1 degrees, where the shape of the ACF is very close to a 
power law. However, it is evident in Figure [5] that the ACF de- 
parts from this model at scales greater than ~ 0. 1 degrees in red- 
shift slices with z p hot > 0.7. The excess of power at large angular 
scale could be explained by the presence of another component 
in the correlation function or of the presence of residual inhomo- 
geneities in the field. We will explore this point more in detail in 
a future paper as it does not impact the results presented here. 
The ACF for the different redshift slices are plotted in Figure [5] 
The best-fitted values of the angular correlation amplitude A w , 
the correlation function slope y and the correlation length ro that 
we derive with the Limber deprojection technique, are reported 
in Table [2] 

Given the uncertainties on the ACF measurements, we adopt a 
fixed slope y = 1.8. This is supported by the fact that when let- 
ting both ro and y vary in the ACF fitting, we find y values very 
similar to 1.8, while the errors on ro increase slightly. 
The values of the correlation length that we derive for the 
3.6pm and 4.5pm photometric samples do not evolve with red- 
shift up to Zphot - 1 .7. In the case of the 3.6pm sub-sample, the ro 
values are in good agreement with the values found in the spec- 
troscopic redshift sample given the size of the error bars. We 
observe slightly larger values of the correlation length for the 
4.5pm photometric redshift sub-sample than for the 4.5pm spec- 
troscopic redshift one. We interpret this as an effect of the in- 
completeness of the spectroscopic samples due to the I-band se- 
lection. To test this we have measured on the photometric sample 
the clustering of the red population which is missed by the spec- 
troscopic sample in each of the 3.6pm and 4.5/im bands. As most 
of the missed population in both cases is at redshifts z P hot > L 
we have been able to perform this analysis only in the redshift 
range Zphot - [1, 1.7] where the available number of galaxies is 
large enough for a clustering signal to be measured. Results are 
presented in Figure [3] We find that the clustering of the missed 
red population in the 3.6 microns spectroscopic sample has a cor- 



T 



a VVDS-Deep (Le Fevre et al. 2005a) 

■ VVDS-SWIRE 3.6 microns spec. 

□ VVDS-SWIRE 3.6 microns phot. (7=1.8) 

• VVDS-SWIRE 4.5 microns spec. 

O VVDS-SWIRE 4.5 microns phot. (7=1.8) 




0.2 0.4 0.6 0.; 



1.2 1.4 



Fig. 6. Clustering evolution of VVDS-SWIRE samples of near- 
infrared galaxies compared to the VVDS-Deep I-band measure- 
ments. Filled symbols refer to measurements from spectroscopic 
samples while open symbols refer to measurements from photo- 
metric samples. 



relation length of ro = 5.6 + 1 .0 only marginally larger than that 
of the bulk of the population, and therefore the clustering sig- 
nal of the photometric and spectroscopic samples are found to 
be very similar. For the red population in the 4.5 microns sam- 
ple, we find that ro = 8.8 ± 2.4, a significantly larger value that 
the main population. This explains why the spectroscopic sam- 
ple slightly underestimates the clustering signal when compared 
to the photometric sample. 

4. Discussion and conclusions 

We first compare our measurements with previously reported 
results on samples of g alaxies selected using Spitzer-IRAC 
bands. lOliveretalJ (12004) use a similar type of galaxy selection, 
as they define a flux-limited sample of SWIRE galaxies with 
S 3.6 > 32pJy, and they compute the ACF and deduce ro using 
Limber deprojection assuming a parametrized redshift distribu- 
tion. The selection criterion they consider is different from the 
one we use, i.e. S3.6 > 9.2pJy. This selects brighter galaxies 
than i n our sample, that are expected to be more strongly clus- 
tered (IPollo et alj2006tlCoil etail2006l) . They find a correlation 
length r = 4.4 + 0.1 h~ l Mpc at z 0.75 with y = 1.8, which 
is higher than our value of ro = 3.7 + 0.4 hr l Mpc at z — 0.81. 
Taking into account that galaxies in their sample are 1.4 mag- 
nitudes brighter and that brighter galaxies are on average more 
clustered, the two results are found to be in good agreement. 
We then compare our results with the clust ering measurements 
of I-band selected VVDS-Deep galaxies (Le Fevre et al. 2005a) 
in the same area. This comparison is shown in Figure [6] where 
the correlation length for the different samples is plotted. We 
observe that the correlation length is globally higher for near- 
infrared selected galaxies compared to optically-selected ones 
and that it is roughly constant over the redshift range z = 
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Table 2. VVDS-SWIRE photometric redshift sample: sub-sample properties and associated measurements of the ACF amplitude, 
the slope and the correlation length. 



Photometric redshift sample 
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redshift interval 


Mean redshift Number of g 


alaxies A w x 10 j 


1 


r (h 1 Mpc) 
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[1.0,1.3] 
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Fig. 5. The angular correlation function as a function of redshift, for the 3.6/j.m and 4.5/im photometric redshift sub-samples. The 
solid lines show the best-fit model when letting the two parameters free, while the dashed lines indicate the best-fit model with fixed 
7 = 1.8. 



[0.2,2.1]. By comparison, ro from the I-selected sample shows 
a steady increase with redshift. As discussed in Section 3, the 
clustering of the 4.5/im-selected galaxy population is higher as 
measured from the photometric sample than from the spectro- 
scopic sample because the latter misses a significant fraction of 
red galaxies due to the additional 17.5 < Iab < 24 cut applied 
by the VVDS-Deep survey selection. When comparing the mea- 
surements at 3.6pnn and 4.5 pun obtained from the photometric 
sample (open symbols in Figure 6), we observe that the cluster- 
ing of the sample selected at the reddest wavelength is slightly 
stronger, although comparable given the errors. This is expected, 



as selecting from longer wavelengths we are sampling preferen- 
tially more early-types than later ones as described below. The 
color dependency of the clustering of galaxies in the VVDS- 
SWIRE sample will be reported elsewhere (de la Torre et al., in 
preparation). 

In order to understand the systematically higher values of clus- 
tering that we find compared to those from optically-selected 
galaxies and to quantify the type of galaxy population included 
in our near-infrared samples, we derive the spectral types of our 
spectroscopic redshift sample galaxies us ing the galaxy classi- 
fication introduced by IZucca et al] d2006l) . This classification is 
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VVDS l m selected VVDS-SWIRE 3.6/j.m VVDS-SWIRE 4.5/im 



VVDS-Deep (Le Fevre et al. 2005a) 
VVDS-SWIRE 3.6 microns spec. 

VVDS-SWIRE 4.5 microns spec. 





0.2 0.4 0.6 0.. 



1 1.2 1.4 1.6 1. 



0.5 1 1.5 0.5 1 1.5 0.5 1 1.5 



Fig. 7. Trend of the linear bias as a function of redshift for the 
VVDS-SWIRE 3.6/im and 4.5/im spectroscopic redshift samples 
compared to the VVDS-Deep I-band one. The dashed curves 
show the theoritical linear halo bias evolution for dark matter 
halos of mass greater than 10 10 , 10 11 , 10 12 , 10 125 and 1O 13 /i"'M 
(from bott om to top). The linear g alaxy bias is computed as de- 
scribed in Mag liocchetti et ail d2000t) from the ro and y measure- 
ments. The theo retical curves have been computed using the fit- 
ting functions of Sh eth & Tormenld 19991) for the halo mass func- 
tion and bias. 



based on the match to an e mpirical set of spec tral energy dis- 
tributions (as described in lArnouts et alj IT999) of the galaxy 
rest-frame colors, obtained from the multi-wavelength informa- 
tion and spectroscopic redshift. The redshift distributions per 
galaxy types are presented in Figure[8]for the I-band, 3.6/-im and 
4.5/im selected samples. The four spectral types defined as type 
1 to 4 correspond respectively to E/S0, early spiral, late spiral 
and irregular/starburst galaxies. One can see in Figure [8] that 
the VVDS-Deep I-selected sample is dominated by late types 
(type 3 and 4) whereas the 4.5/im selected sample is dominated 
by early types (type 1 and 2). In the 3.6/im sample the dif- 
ferent types balance out. Early type galaxies ar e known to be 
more clustered than late type ones at z = (e.g. Norber g et al 



2002: lZehavi et al.l2 005) and also at high redshift dMeneux et a! 



2006), thus the proportion of early type galaxies that we found 
in the different samples can explain why globally we measure 
higher values of clustering when galaxies are selected from red- 
der observed wavelengths. 

On the other hand, understanding the relatively constant cluster- 
ing as a function of redshift that we observe in our near-infrared 
samples is not straightforward. Selecting galaxies at 3.6//m or 
4.5/im, we are probing the rest-frame near-infrared of galax- 
ies, more sensitive to the emission from old stars, and therefore 
our near-infrared selection is mostly driven by stellar mass. As 
dark matter halo clustering is expe cted to decrease with increas- 
ing redshift dWeinberg et al.| [2004). and under a reasonable as- 
sumption that stellar mass traces the underlying dark matter halo 
mass, one would then expect the correlation length of our sam- 



Fig. 8. Redshift distribution per spectral types for the 3.6/im and 
4.5/im spectroscopic redshift sub-samples compared to the 
VVDS-Deep samp le. We use the spectral types defined in 
IZucca et alj 12006): type 1 corresponds to E/S0, type 2 to early 
spiral, type 3 to late spiral and type 4 to irregular/starburst galax- 



ies. 



pie to decrease with increasing redshift. On the contrary, as we 
select intrinsically more luminous galaxies with increasing red- 
shift because of the magnitude selection of the sample, and since 
more luminous galaxies are more clustered dPollo et al.1 12006; 
ICoif et al.l2006l) . we expect higher clustering strength as redshift 
increases. We suggest that these two competing effect, working 
simultaneously, combine to produce the roughly constant clus- 
tering amplitude with redshift that we observe. 
In order to check this hypothesis we look at the evolution of 
the linear galaxy bias for the differe nt samples. We compute the 
galaxy linear bias bi as described in Magli occhetti et al"] ([2000) 
from our ro and y measurements. We define b^ as, 



bdz) = 



A /c y (r (z)/8) y 

CT &,mass(z) 



72 



(3-y)(4-y)(6-y)2? 



(8) 



This assumes a linear evolution of the rms mass fluctuations with 
cosmic time, i.e. o-^ mass (z) = o-^ mass (z = 0)D(z), D(z) being the 
linear growth factor. We compare it to the theoretical halo bias 
evolution for dark matter halos of constant mas s. To evaluate the 
linear halo bias we use the fitting functions of I Sheth & Tormenl 
(1999) for the halo mass function and bias. As shown in Figure 
Qthe linear galaxy bias evolves more rapidly than the linear halo 
bias, indicating that these galaxies are not strictly mass-selected. 
Indeed because of our magnitude selection, we are probing more 
luminous, hence more biased galaxies at increasing redshifts. 
However, at least up to z - 1 .5 (and given the uncertainties), the 
linear galaxy bias for the rest-frame near-infrared galaxy sam- 
ples seems to stay in between the linear halo bias tracks for halos 
of mass greater than 10 12 and 10 12 5 h~ l M G . This spread in mass 
is lower than the one for the rest-frame UV-optical VVDS-Deep 
sample, making this sample more suitable to trace the passive 
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growth of clustering in the mass. This comparison of the linear 
galaxy bias and the linear halo bias is consistent with the idea 
that the clustering dependence on luminosity and stellar mass 
combine to produce a relatively constant clustering with redshift 
as suggested in the previous paragraph. 

It is particularly interesting to note that at z — 1 -5, the correlation 
length of galaxies selected from the UV-optical rest-frame and 
hence actively star-forming is very similar to that of galaxies se- 
lected in the rest-frame near-infrared dominated by older stellar 
populations. This is quite different from what is observed in the 
local Universe where selecting gal axies in the K s band from the 
2MASS survey (iMaller e t al. 2005]), it is found that these galax- 
ies are significantly more clu stered than those sele cted at opti- 
cal wavelengths (SDSS EDR, IConnollv et al.ll2002b . The galax- 
ies which are making stars at z — 1.5 therefore reside in sim- 
ilarly clustered regions as those containing the bulk of stellar 
mass, while at z — the star formation happens in a population 
which is much less clustered than that where the bulk of stel- 
lar mass already resides. This result is consistent with the ob- 
served flattening of t he color-density relation above z - 1.2-1.5 
(Cucciati et al. 2006), which indicates that at these redshifts, 
blue star forming galaxies and red older galaxies are found to 
reside with equal probability in high or low density regions. 
Furthermore, we note that the clustering of massive luminous 
galaxies at z — 1.5 is measured to be rrj 4 lr l Mpc, indepen- 
dently of the selection at optical or near-infrared wavelengths, 
while at lower redshifts, the difference in the clustering strength 
measured from these samples increases. Galaxies selected from 
their UV-optical rest-frame show lower and lower clustering, 
whereas galaxies selected in the near-infrared (i.e. more mass- 
selected) stay equally clustered. Taken together, this can be in- 
terpreted as evidence for star formation shifting from high den- 
sity to low density regions, as cosmic time increases, another 
manifestation of the downsizing trend whic h has been indicate d 
by observations over the last few years (e.g. lCowie et al. I fl996h . 
This underlines the pri mary role of the ma ss in regulating star 
formation in galaxies (Gava zzi et al.l 12002b : while at high red- 
shift star formation was strong in massive and more clustered 
galaxies, today it is mainly limited to lower mass galaxies more 
uniformly distributed. 
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